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The search for Majorana bound states in solid-state physics has been limited to materials which 
display a gap in their bulk spectrum. We show that such unpaired states appear in certain quasi- 
one-dimensional Josephson junctions arrays with gapless bulk excitations. The bulk modes mediate 
a coupling between Majorana bound states via the Ruderman-Kittel-Yosida-Kasuya mechanism. 
As a consequence, the lowest energy doublet acquires a finite energy difference. For realistic set of 
parameters this energy splitting remains much smaller than the energy of the bulk eigenstates even 
for short chains of length L ~ 10. 


Introduction. An intensive search for Majorana 
fermions [T] is underway in solid-state devices [2]. The 
vast majority of the proposals consist in zero energy 
boundary modes in materials with a gaped bulk spec¬ 
trum. A spin-less superconducting wire or topological 
insulator in two or three dimensions fall in this cate¬ 
gory [3H5]. 

We propose an alternative approach for the observa¬ 
tion of Majorana fermions in Josephson Junctions Arrays 
(JJA). We will show that certain quasi-one-dimensional 
JJA can display Majorana zero modes at their bound¬ 
aries. These modes are protected from mixing with 
higher energy excitations although bulk spectrum is not 
gapped. The existence of low energy Majorana could 
then be proved by spectroscopy lam- 

in this letter, we first explain the JJA system and how 
to model it with an Ising-like Hamiltonian. Then, a qual¬ 
itative argument is employed to obtain the low-energy 
effective theory using Majorana boundary modes. Nu¬ 
merical results will confirm the validity of this effective 
theory and show that Majorana modes are indeed pro¬ 
tected. Finally, we discuss problems that may arise in 
the experimental realization of our proposal. 

Experimental set-up We consider three identical lad¬ 
ders of Josephson junctions coupled together as in Fig. 
&) i- Each ladder has a unit cell with ’Targe” 
and ’’small” junctions. Fig. &)• Their corresponding 
Josephson energies are Ejl and Ejs, where Ejl > Ejs. 
The three ladders are closed at the ends by a junction 
with Josephson energy Eje. We assume that charging 
energies are much smaller than Josephson ones for all the 
junctions. All the closed circuit in the ladders are at full 
frustration, that is, they are threaded by a magnetic flux 
equal to half of the flux quantum. The two larger loops 
are threaded by magnetic fluxes ipi and ip 2 , respectively. 
For ladder p, the phase difference in the left vertical and 
horizontal red junctions in the loop k are denoted by 
4>p{fJk) and (j)p{2k -|- 1), respectively (see Fig.[^) . 

At full frustration, the properties of one infinite chain 
are invariant under translation by one small closed loop 



FIG. 1. (a) Three ladders of Josephson junctions coupled 

together. The two large loops are threaded by magnetic fluxes 
i/Ji and ifi 2 - (b) Unit cell of one ladder. Josephson energy 
of ’’large” junction (red) is larger than the one of a ’’small” 
junction (green), Fjl > Ejs. All the cells are threaded by 
a magnetic flux equal to half of the flux quantum, (c) Loop 
at the left boundary (the one at right is symmetric). The 
Josephson coupling for the blue junction at the chains end is 
denoted as Eje ■ (d) For small temperature and near a phase 
transition, the three Josephson ladders system maps to three 
Ising chains with transverse held (solid lines) and coupling 
between their ends (dotted lines). Majorana zero modes are 
located at the chains boundaries (dark regions). 


and reflection through horizontal axes. We analyze the 
energy of one unit cell as a function of the phase differ- 
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ence in the ’’large” junctions: 

Up^k = -Ejl I cos [^f)p(2fc)] + cos [(j)p{2k + 1)] | 

EEjg cos \(j>p{2k) + cj)p(2k + 1) + (j)p(2k + 2)], (1) 

where Up^k is the energy of the loop k of ladder p. The 
total energy has a global Z 2 symmetry given by = 
—4)p{i) in each junction i. In the regime Eji, ^ Ejs, 
the ground-state corresponds to all phases equal to zero 
4>p{i) — 0 and the symmetry is preserved. However, a 
broken symmetry phase occurs for Ej^ ^ Ejg, as in 
this regime the ground-state acquires a ^p(j) ^ 0. The 
critical point is located near Ejl/Ejs ~ 5 [S]. 

The specific details of the system are irrelevant near 
the phase transition and the properties of one ladder are 
described by an Ising chain with transverse field. The 
mapping can be written explicitly by taking each ’’large 
junction” as a 1/2-spin with component ^1 — 4>p{i) in 
the direction of the field. Near the phase transition, 
only the lowest non-zero order in the phase differences 
in ’’large” junctions are relevant. Then, the first term in 
Eq. 0 represents the field contribution and the second 
the Ising coupling between nearby spins. 

The junctions located at the ends of the chains, blue 
crosses in Fig. §i). couple the three ladders: 

klc = Eje I cos [01(1) -I- 02(1) + 03(1)] 

-I-cos [01 (L)-I-02 (T)-I-03 (T)] I, (2) 

where L is the double of the number of squares in each 
ladder. This contribution gives an extra coupling be¬ 
tween spins at the boundaries in the Ising model. 

Effective model. Near the phase transition, the 
Hamiltonian of the three ladders of JJA maps to three 
Ising chains with transverse magnetic field and coupling 
between their ends H = Hp + Elc, where: 

L-l L 

Hp=-jY,<^;ikK{k + l)-hY,<^pik), ( 3 ) 

k^l k^l 

H,=- (4) 

P=1 

The index p = 1, 2,3 is the chain number and it has peri¬ 
odic boundary conditions. The length of each chain is L 
and thus, the total number of sites of the three chains is 
N = 3L. The operators are Pauli ma¬ 

trices acting on the Hilbert space of the spin at site i of 
chain p. The parameters J and h control the usual Ising 
interaction and magnetic field in the transverse axis, re¬ 
spectively. They can be related to the original Josephson 
coupling energies as J/h ~ Eje/Ejs The value of Jc sets 
the strength of the coupling between different chains. It 
can be controlled by the junctions located at ends of the 


Josephson ladders Jc/J ^ Eje/Eje- We are interested 
in the regime /i/J ~ 1 and Jc/J ~ 1- A schematic rep¬ 
resentation of the Hamiltonian appears in Fig. ©)• 

In Eqs. ( |I|2[ ), we have implicitly assumed that quan¬ 
tum fluctuations in the phase of the ’’small” junctions are 
zero. In practice, these fluctuations are small but non¬ 
zero. We have also neglected any noise in the flux thread¬ 
ing superconductor loops. These two contribution repre¬ 
sent sources of incoherent noise for our effective Hamilto¬ 
nian. We will comment later on how those contributions 
affect our results. 

We express the low-energy degrees of freedom for 
Hamiltonian Eq. in terms of Majorana fermions. 

We use a multichannel version of the Jordan-Wigner 
transformation [3] . This mapping requires to enlarge the 
Hilbert space with an extra I/2-spin. Operators acting 
on this spin are denoted by cr’“(0), cr^(O) and cr^(O). The 
original spin operators are mapped to fermions via: 


cUk) = Pp{-l)^i«^' 


{k ), 


( 5 ) 


where np{j) = [o'|(p) -I- l] /2. The rjp acts on the added 
1/2 fermion. They are: 

771 = a"(0)(-l)^^+^^ 

V2 = <jH0) 


where Np is the total number of fermions in chain p. The 
operators defined in Eq (§ fulfill the fermionic algebra. 
Finally, Majorana fermions are defined as tpp{2k — 1) = 
Cp{k) + cj,(k) and il}p(/2k) = i \cp{k) — cj,(/c)]. 

We focus on the regime in which each of the Ising 
chains are critical, J = h. In this case, our model can be 
written as: 

2L 

iH = J E ^(j)0(j +1)+ (6) 

1=1 

Jc^eabcVbVc [06(1)00(1) -f 06(2T)0 c(2L)(-1)^'’+^"] , 

abc 


where = (0i(/), 02(/)_j 03(j))^ The coupling at 

the left boundary is S ■ (0(1) x 0(1)), where 5*“ = 
{i/2)eabcVaVb- If only this coupling is considered, the 
Hamiltonian describes an over-screened two channel 
Hondo (2CK) model (TUHIS]- The role of the impurity 
is played by the 1/2-spin S and the relevant degrees of 
freedom of the conduction electrons are described by the 
bulk spins of the three Ising chains. At temperatures be¬ 
low the Hondo temperature Tk ~ J exp (—const J/Jc) the 
two-channel Hondo model displays the universal quan¬ 
tum critical behavior. At the Quantum Critical point 
(QCP) of the 2CH model possesses a finite ground state 
entropy In \/2 which originates from the presence of a 
zero energy Majorana mode. This mode presents what is 
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left from the impurity spin S=l/2. The leading irrelevant 
operator at the QCP describes a coupling of this mode 
to the product of three bulk Majorana at the impurity 
point. 

It stands to reason that at L —>■ oo model (§ describes 
two independent 2CK models. Formally such decoupling 
is achieved by declaring (—l)'^“77a as an independent Ma¬ 
jorana fermion. Then its anticommutator with rja is zero 
on average. At finite L Hamiltonian ®> possesses an ad¬ 
ditional energy scale generated by a tunneling between 
the Majorana modes located at different ends of the sys¬ 
tem. Below this scale one has to see deviations from the 
two-channel Kondo physics. Although model § is not 
integrable, we can analyze qualitatively the low energy 
theory starting from the strong coupling limit Jc/J — 0. 
There are Majorana fermions at each end of the chains, 
Ml = '0i(l)^2(l)'*/'3(l) and Mit = '0i(2T)^/’2(2T)^/’3(2T), 
which commute with the strong coupling Hamiltonian. 
That implies a two-fold degeneracy of the whole spec¬ 
trum. Indeed, the occupancy of the fermionic state 
/ = [ml — does not affect the energy. We notice 
that the fermion / is constructed from Majorana at dif¬ 
ferent ends of the chain. 

If one adopts the hypothesis of independence of two 
2CK QCPs at L —)■ oo and consider the finite size effects 
as a perturbation their contribution to the low energy dy¬ 
namics can be extracted from the effective Hamiltonian 
computed as m- 

Hr r-^/"[/ZLV'l(2)V'2(2)V'3(2) + 

MflV'i(2A-l)V'2(2L-l)V'3(2T-l)]. (7) 


We expect that Majorana at each end couple via RKKY 
mechanism through the low energy excitations [HHiT]. 
Indeed, integrating over the fermions in 0 we obtain the 
effective Hamiltonian for the Majorana modes: Heff = 
t where T ^ f G{L, t)^dt ^ T^^{vf/LY■ The 
Green function of the bulk Majorana is 


G{x,t) =e 


„ikpx 




2L sin 


+e 


— ikpx 


7r(x + itvF)/‘2L 
vf'x 


2Lsin 


7r(a; — itvF)/2,L 


( 8 ) 


where = 7r/2 is the Fermi wave vector and vf is the 
Fermi velocity m- In this simpler case, the degener¬ 
acy of the spectrum is lifted by a factor T ~ L~^. The 
effective Hamiltonian can be written as 


if f--). 


( 9 ) 


If this result holds for the full model Eq. (3j4), the ex¬ 
istence of Majorana fermions in the JJA discussed above 
would be signaled by a mode with energy going as L~^. 
The bulk eigenstates correspond to that of the Ising chain 
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FIG. 2. Difference between the energy of the first excited and 
ground states, Ei — Eo, of the three Ising chains Eqs. ( |3|4| | as 
a function of the inverse of the chains length 1/L. Each chain 
is at the critical point h = J. The strength of the coupling 
between chains is Jc = J. The size of the system goes from 
N — 6 to N = 21, where N = 3L. The line fits the data 
to ax^ + b, where a, b and p are free parameters. The result 
gives an exponent p — 2.05 ± 0.05. The other parameters of 
the fit are b = (1.0±0.5) x 10"'‘ and a = (158.7±0.6) x 10"^ 
Inset: the black points are Ei — Eq as a function of 1/L, 
as in the main panel. The red squares correspond to the 
difference between energy of the third excited and ground 
states, E 3 — Eq. 


with transverse field at the critical point |18j . whose en¬ 
ergy levels scale as e > 1/L. Then, for sufficiently long 
chains the unpaired Majorana are protected. We are go¬ 
ing to check that this is actually the case for the full 
model Eqs. (3j4| with parameters J = h = Jc- In fact, it 
turns out that for the parameters chosen the low energy 
doublet is still protected even in the case of short chains. 


Numerical results. We employ numerical methods to 
show that our previous result based on Majorana bound¬ 
ary modes is also valid for the full model Eqs. 
Specifically, exact diagonalizations of the three Ising cou¬ 
pled chains are carried out for sizes ranging from N = 6 
to N = 24. We use ARPACK libraries as only a few 
eigenstates are required [H]. These libraries are based 
on Lanczos algorithm which is specially suited to per¬ 
form efficient partial diagonalization of sparse matrices. 


In Fig. the difference in energy between the first 
excited and ground-state, Ei — Eq, appears as a function 
of 1/L. The data have been fitted to a function a*xP + b, 
where a,b,p are free parameters. The result is the solid 
black line with p = 2.05±0.05 and b = (1.0±0.5) x 10“"^. 
These values are in excellent agreement with our previous 
analysis based on RKKY interactions between Majorana 
particles. 


The low-energy doublet can be resolved in experiments 
only if higher excited states have much larger energy. Of 
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FIG. 3. Difference between the energy of the first excited and 
ground states, Ei — Eq, as a. function of the inverse of the 
chains length 1/L. A random magnetic field in the x-axis, 
following a uniform distribution in the interval [—IF/2, IF/2], 
is added to the original Hamiltonian with h = J = Jc- Each 
set of data corresponds to a different value of IF, as indicated 
in the legend. The lines fit the points for each disorder to a 
function ax^ + b, where a, b and p are free parameters. The 
value of the fitting at the origin range approximately from 
b = 0.014 to fo = 0.0004 and the exponent from p = 2.6 to 
p = 2.07. 


course, this occurs in the thermodynamic limit, as the 
doublet energy Ei — Eq ^ 1 jL?' is much smaller than the 
energy of the second excited state E^ — E^ ^ XjL. The 
red squares in the inset of Fig. [^represent the difference 
in energy between the second-excited and ground-state, 
E 2 — Eq, as a function of 1/L. The black dots are Ei — Eq 
as in the main panel. We can see that E 2 — Eq is more 
than two orders of magnitude larger than Ei — Eq. The 
size of the chains that we have employed runs from L = 2 
to L = 8. We conclude that Majorana modes are well 
protected even for small chains. 

We discuss the source of incoherent dynamic in the 
experimental set-up. Fluctuations due to the quantum 
nature of the junctions can be kept small by employing 
a large ratio Ej/Ec > 100 [8]. In this case, phase slips 
are neither important as their probability decrease expo¬ 
nentially in Ej/Ec [7]. 

The largest source of noise is caused by time-dependent 
fluctuations in the magnetic field threading the large 
loops (fi, ip 2 (Fig. Si). which produce additional con¬ 
tributions in the phases of ” small” junctions of the order 
of ~ ^i/L. These can be taken into accounts by intro¬ 
ducing a small random argument in the third cosine of 
Eq. 0 or, in the language of our effective spin model, by 
adding a small random magnetic field in the x-axis. It is 
important to realize that this noise can be minimized by 
constructing the JJA with a small area between nearby 
chains. We are going to show that Majorana modes are 


still protected if this noise is kept small enough. 

We present results for exact diagonalization of our 
model in the presence of a random magnetic field in 
the x-axis. In Fig. the energy difference of the first 
two levels ALi is computed as a function of 1/L for the 
three Ising chains, Eq. (3]4), plus a small perturbation 
ATJa; = where ep(fc) are randomly dis¬ 

tributed in the interval [—IF/2, IF/2]. The solid lines fit 
each set of data to a function a * + 6, with free pa¬ 

rameters a,b,c. For the largest disorder IF = 2 * 10“^ J, 
the fitting function provides p = 2.6 and it does not go 
trough the origin b = 0.014. This is expected as the typ¬ 
ical energy of the perturbation, IF, is roughly the same 
as the one of the first level AE for the largest size. The 
points for the smaller disorder display a more similar be¬ 
havior to the one in the clean system. Furthermore, the 
low energy mode remains in all cases well separated from 
the bulk spectrum in the presence of small disorder. 

Summary. We have shown that Majorana modes ap¬ 
pear in a system composed by three ladders of Joseph- 
son junctions. The signature of Majorana is a low en¬ 
ergy mode with energy going as ^ 1/L^, while the bulk 
spectrum eigenstates have energies e ~ 1/L. The lowest 
energy doublet is then protected in the thermodynamic 
limit. We have further checked that the energy difference 
of the lowest doublet is more than two order of magni¬ 
tude smaller than energy of bulk eigenstates for chains of 
length L ^ 10. So, boundary modes are also protected for 
system sizes which are accessible to experiments. Finally, 
we have analyzed the sources of incoherence and noise in 
the proposed experimental set-up. We conclude that our 
JJA can be a good experimental system to search for the 
Majorana particle. 
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